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Abstract 


Large  DND  projects  and  acquisitions  are  exposed  to  more  than  one  foreign  currency  at 
the  same  time  which  complicates  management’s  foreign  exchange  risk  assessments.  We 
extend  the  Centre  for  Operational  Research  and  Analysis’  (CORA)  in-house  Generalized 
Autoregressive  Conditional  Heteroskedasticity  (GARCH)  models  to  a  full  multivariate  set¬ 
ting.  Our  extensions  involve  two  models  types:  multivariate  GARCH  and  copula-GARCH. 
We  find  that  both  models  give  qualitatively  similar  value-at-risk  (VaR)  estimates,  and  that 
both  models  provide  a  much  improved  risk  assessment  relative  to  the  current  practice  -  cor¬ 
recting  VaR  estimates  on  the  order  of  25%  in  cases  in  which  multiple  currency  exposures 
are  of  similar  size.  Using  the  USDCAD,  the  EURCAD,  and  the  GBPCAD,  we  demon¬ 
strate  estimation  techniques  for  each  model.  Finally,  we  show  the  strength  of  our  improved 
models  through  a  100-day  VaR  calculation. 
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Significance  for  defence  and  security 


The  Department  of  National  Defence  (DND)  expends  approximately  10%  of  its  annual 
budget  in  foreign  currencies  in  support  of  equipment,  acquisitions,  and  operations.  The  for¬ 
eign  currency  transactions  within  DND’s  budget  expose  the  department  to  financial  market 
risk  in  the  form  of  fluctuating  exchange  rates,  placing  additional  risk  mitigation  pressures 
on  internal  management.  Given  the  Government  of  Canada’s  planned  military  acquisitions 
in  the  form  of  the  CF-188  replacement  and  the  renewal  of  the  Royal  Canadian  Navy  -  both 
of  which  will  contain  significant  foreign  content  -  exchange  rate  risk  has  quickly  become 
an  important  risk  element  for  DND  management  to  understand. 

In  this  study,  we  extend  our  in-house  financial  risk  models  to  help  DND  understand  the 
extent  of  risk  mitigation  that  occurs  naturally  through  exposure  to  more  than  one  currency. 
DND’s  main  foreign  currency  exposure  involves  the  US  dollar,  the  euro,  and  the  British 
pound  sterling.  To  gauge  the  effect  of  adverse  currency  fluctuations  on  DND’s  budgets  and 
planning,  we  require  an  understanding  of  the  strength  of  natural  diversification  benefits 
within  value-at-risk  (VaR)  analyses.  We  examine  two  methods  to  attack  this  problem.  Our 
first  approach  uses  multivariate  correlation  methods  which  extract  time  dependence  and 
persistence  in  correlations  between  exchange  rate  pairs.  In  our  second  approach,  we  build 
dependency  models  between  exchange  rate  pairs  using  copulas.  Both  modelling  techniques 
offer  decision  makers  a  deeper  understanding  of  portfolio  VaR  -  correcting  the  indepen¬ 
dence  assumed  VaR  estimate  by  25%  in  cases  in  which  multiple  currency  exposures  are  of 
similar  size. 
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Resume 


Dans  ses  grands  projets  et  achats,  le  MDN  doit  souvent  composer  avec  plusieurs  devises 
a  la  fois,  ce  qui  complique  son  evaluation  des  risques.  Nous  avons  amplifie  des  modeles 
d’heteroscedasticite  conditionnelle  autoregressive  generalisee  (GARCH)  du  Centre  d’ ana¬ 
lyse  et  de  recherche  operationnelle  (CARO)  de  maniere  a  y  inclure  toutes  les  variables 
necessaires.  Nous  avons,  de  fait,  modifie  deux  types  de  modele  :  le  modele  GARCH  a 
variables  multiples  et  le  modele  GARCH  a  copules.  Nous  avons  constate  que  les  deux  mo¬ 
deles  donnent  des  estimations  de  la  valeur  a  risque  (VaR)  de  qualites  comparables  et  qu’ils 
permettent  d’effectuer  une  meilleure  evaluation  du  risque  qu’avec  la  methode  courante.  On 
obtient  desormais  des  estimations  de  la  VaR  25%  plus  exactes  pour  les  cas  ou  plusieurs 
devises  entrent  en  jeu.  Au  moyen  des  combinaisons  USD-CAD,  EUR-CAD  et  GBP-CAD, 
nous  avons  demontre  les  techniques  d’estimation  des  deux  modeles.  Enfin,  nous  avons  de- 
montre  la  fiabilite  de  nos  modeles  ameliores  au  moyen  d’un  calcul  de  la  VaR  sur  100  jours. 
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Importance  pour  la  defense  et  la  securite 


Le  ministere  de  la  Defense  nationale  (MDN)  depense  environ  10%  de  son  budget  annuel 
en  devises  etrangeres  pour  l’entretien  de  son  equipement,  ses  achats  et  ses  operations.  Ces 
transactions  en  devises  etrangeres  exposent  le  MDN  a  un  risque  financier  lie  aux  varia¬ 
tions  du  taux  de  change,  et  les  responsables  de  la  gestion  interne  se  voient  done  presses 
de  trouver  des  strategies  d’ attenuation.  Comme  le  gouvemement  du  Canada  prevoit  faire 
d’importantes  acquisitions  militaires,  notamment  pour  remplacer  le  CF188  et  renouveler  la 
Marine  royale  canadienne  (deux  projets  dans  lesquels  il  fera  beaucoup  affaire  avec  l’etran- 
ger),  le  risque  lie  au  taux  de  change  est  rapidement  devenu  un  element  important  que  la 
direction  du  MDN  se  doit  de  maitriser. 

Dans  notre  etude,  nous  avons  amplifie  nos  modeles  maison  d’ evaluation  du  risque  finan¬ 
cier  dans  le  but  d’aider  le  MDN  a  comprendre  l’importance  de  P  attenuation  des  risques 
lorsqu’on  doit  composer  avec  plus  d’une  devise.  Les  principals  devises  etrangeres  que  le 
MDN  rencontre  le  plus  frequemment  sont  le  dollar  americain,  l’euro  et  la  livre  sterling 
britannique.  Pour  bien  mesurer  les  effets  negatifs  que  peuvent  avoir  les  fluctuations  mo- 
netaires  sur  le  budget  et  la  planification  du  MDN,  il  faut  connaitre  le  poids  des  benefices 
de  la  diversification  naturelle  dans  les  analyses  de  la  valeur  a  risque  (VaR).  Pour  aborder 
ce  probleme,  nous  avons  examine  deux  methodes.  Dans  la  premiere,  nous  avons  utilise 
des  modeles  de  correlation  a  multiples  variables  qui  isolent  la  dependance  et  la  persistance 
temporelles  dans  les  correlations  entre  deux  devises.  Dans  la  seconde  methode,  nous  avons 
elabore  des  modeles  de  dependance  a  copules.  Ces  deux  methodes  offrent  aux  decideurs 
une  meilleure  comprehension  de  la  VaR  d’un  portefeuille,  ce  qui  leur  permet  d’obtenir  des 
estimations  25%  plus  exactes  dans  les  cas  ou  plusieurs  devises  entrent  en  jeu. 
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1  Introduction 


1.1  Background 

The  Department  of  National  Defence  (DND)  expends  approximately  10%  of  its  annual 
budget  in  foreign  currencies  in  support  of  equipment,  acquisitions,  and  operations.  The  for¬ 
eign  currency  transactions  within  DND’s  budget  expose  the  department  to  financial  market 
risk  in  the  form  of  fluctuating  exchange  rates,  placing  additional  risk  mitigation  pressures 
on  internal  management.  Given  the  Government  of  Canada’s  planned  military  acquisitions 
in  the  form  of  the  CF-188  replacement  and  the  renewal  of  the  Royal  Canadian  Navy  - 
both  of  which  will  contain  significant  foreign  content  -  foreign  exchange  risk  has  quickly 
become  an  important  risk  element  for  DND  management  to  understand. 

Over  the  last  seven  years,  the  Centre  for  Operational  Research  and  Analysis  (CORA)  has 
developed  in-house  financial  risk  management  tools  to  help  DND  understand  foreign  ex¬ 
change  market  risk  within  its  budgets  [1],  [2].  CORA’s  approach  uses  two  main  lines 
of  inquiry  -  models  based  on  Generalized  Autoregressive  Conditional  Heteroskedasticity 
(GARCH)  [1],  and  models  based  on  the  internal  consistency  of  foreign  exchange  derivative 
prices  [2].  Both  modelling  methods  yield  value-at-risk  (VaR)  [3]:  given  a  confidence  level 
p,  and  a  prescribed  time  horizon  T,  the  VaR  yields  the  foreign  exchange  loss  level  that 
is  exceeded  with  probability  1  —  p  at  time  T.  Decision  makers  can  use  VaR  calculations 
matched  to  foreign  payment  schedules  to  assess  the  likelihood  of  project  budget  shortfalls. 

In  conjunction  with  risk  modelling,  CORA  has  helped  DND  prioritize  corporate  risk  by 
using  expert  opinion  with  rank  correlation.  The  analysis  showed  that  “increased  financial 
pressure  due  to  external  factors  represents  one  of  the  highest  risk  factors  among  all  corpo¬ 
rate  activities  at  DND”  [4].  Specifically,  ADM(Fin-CS)  and  ADM(Mat)  singled  out  foreign 
exchange  risk  as  a  significant  liability.  Senior  decision  makers  at  DND  continue  to  hold 
the  view  that  foreign  exchange  represents  a  significant  risk  in  DND  planning  activities. 

In  this  study,  we  extend  our  in-house  GARCH  models  to  help  DND  understand  the  ex¬ 
tent  of  risk  mitigation  that  occurs  naturally  through  exposure  to  more  than  one  currency. 
DND’s  main  foreign  currency  exposure  involves  the  US  dollar,  the  euro,  and  the  British 
pound  sterling.  To  gauge  the  effect  of  adverse  currency  fluctuations  on  DND’s  budgets  and 
planning,  we  require  an  understanding  of  the  strength  of  natural  diversification  benefits 
within  VaR  analyses.  We  examine  two  methods  to  attack  this  problem.  Our  first  approach 
uses  multivariate  GARCH  models  which  extract  time  dependence  and  persistence  in  corre¬ 
lations  between  exchange  rate  pairs.  While  this  method  offers  a  powerful  modelling  tool, 
the  method  implicitly  contains  restrictive  assumptions  about  the  nature  of  the  dependency 
between  the  two  exchange  rates.  We  address  these  shortcomings  in  our  second  approach  by 
building  dependency  models  between  exchange  rate  pairs  using  copulas.  Both  modelling 
techniques  offer  decision  makers  a  deeper  understanding  of  portfolio  VaR  in  which  more 
than  one  currency  affects  total  project  costs. 
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1.2  Scope 

ADM(FinCS)  requires  an  understanding  of  foreign  exchange  risk  using  VaR  calculations 
for  its  budgets.  In  this  paper,  we  propose  a  VaR  analysis  based  two  methods: 

•  multivariate  GARCH  models;  and 

•  copula-GARCH  models. 

We  apply  these  techniques  to  EURCAD-USDCAD,  and  GBPCAD-USDCAD  exchange 
rate  pairings  using  daily  data  from  April  30,  2004  to  May  1,  2014. 


2 


DRDC-RDDC-2014-R62 


2  Generalized  autoregressive  conditional 
heteroskedasticity  (GARCH)  basics 


Foreign  exchange  markets  show  no  evidence  of  autocorrelation  in  the  return  data,  implying 
that  past  observations  of  foreign  exchange  prices  do  not  help  forecast  future  returns.  This 
observation  -  that  past  prices  offer  no  predictive  power  and  thus  no  advantage  in  trading  - 
conforms  with  the  efficient  market  hypothesis  (EMH)[5].  However,  even  though  uncondi¬ 
tionally  return  variances  are  nearly  stationary1,  for  most  assets,  and  especially  for  foreign 
exchange,  they  are  not  conditionally  stationary.  To  understand  the  return  sequence  of  for¬ 
eign  exchange  data,  we  must  account  for  time  changing  conditional  return  variances.  The 
framework  for  modelling  this  phenomenon  is  called  generalized  autoregressive  conditional 
heteroskedasticity  (GARCH)  [7]  (see  [8]  for  a  complete  review). 


Under  the  simplest  GARCH  model,  called  GARCH(1,1),  the  variance  process  takes  the 
form, 

crf+i  =  co  +  aRj  +  /3  a f2,  (1) 

where  a  +  /3  <  1  and  Rt  denotes  the  asset  return.  In  this  paper,  we  use  the  arithmetic  return, 


Rt 


5f+ 1  —  St 


St 


5 


(2) 


where  St  is  the  observed  exchange  rate  at  market  close.2  In  GARCH  modelling,  we  assume 
that  the  return  sequence  follows, 

Rt  =  otZt,  (3) 

where  zt  follows  an  iid  process. 


In  addition  to  modelling  univariate  time  series  with  GARCH  models,  we  can  model  the 
linear  correlation  between  time  series  using  similar  methods.  Note  that  we  can  re-write 
the  univariate  GARCH(1,1)  model  by  recognizing  that  the  unconditional  variance  has  the 
relationship, 


o 


2 


(7 


2 


E(ct2+1)  =  co  +  aE(Rj)  +  /3E(of) 

co  +  a  a2  +  /3  a2,  implying, 

CO 

1  -a-/r 


(4) 


where  we  assume  that  the  daily  return  has  a  vanishing  expectation.  Using  this  approach, 
we  can  eliminate  co  from  the  equations,  allowing  us  to  write  the  bivariate  GARCH(1,1) 

'For  details  on  long  term  predictability  in  asset  returns,  see  [6]. 

2In  the  literature,  log  differences  are  often  taken  as  the  return,  but  we  conform  to  the  approach  taken  in 
[9].  Note  that  for  small  changes,  log  and  arithmetic  returns  are  nearly  identical. 
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process  as, 


Qt+ 1  =E(zfz}r)(l  -  a  -  p)  +  atztz?)  +  pQt 

qn,t+i  qn,t+i  A 
qn,t+i  qn,t+i  J 

=(p  0(i-“-/5) 

+  a(  q’f  Zl’f2^+p(qn^ 

V  ZMZ2,t  25, f  J  V  912, r 


912, f 
<722,/ 


(5) 


where 


P/+1  = 


<7l2,f+l 


(6) 


a/<71  l,r+l<?22,r+l 

and  denotes  the  i-th  univariate  GARCH  filtered  return.  The  unconditional  correlation 
coefficient  is  given  by, 


1 


P  =  —  ^2l,/22,i 


(7) 


t=l 


The  multivariate  GARCH  model  allows  us  to  track  the  time  dependence  and  the  persistence 
of  the  correlation  coefficient.  Again,  we  use  maximum  likelihood  methods  for  parameter 
estimation.  The  underlying  structure  of  the  multivariate  GARCH  model  is  a  bivariate  nor¬ 
mal  distribution  with  a  time  changing  correlation  coefficient. 
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3  A  brief  introduction  to  copulas 


In  statistical  modelling,  we  often  desire  an  understanding  of  the  full  joint  distribution  func¬ 
tion  of  a  collection  of  random  variables.  For  example,  the  2-dimensional  case  with  random 
variables,  X  and  Y,  have  the  joint  distribution  function, 

P(*<x,F<y)=//(x,y),  (8) 

and  in  an  optimal  decision  making  context,  we  might  require  an  inference  of  H(x,y)  from 
the  empirical  data.  When  confronted  with  the  problem  of  estimating  the  dependence  be¬ 
tween  random  variables,  we  often  rely  on  key  statistical  parameters,  such  as  the  covariance 
or  the  correlation  coefficient.  Unfortunately,  dependence  runs  far  deeper  than  simple  co- 
variances. 


To  gain  an  appreciation  of  the  problem  at  hand,  suppose  that  we  make  independent  draws 
from  the  standard  normal,  keeping  the  result  and  its  square  -  that  is,  we  construct  the 
bivariate  process  (X .X2).  Clearly  X2  depends  perfectly  on  X,  knowledge  of  the  draw 
implies  complete  knowledge  of  its  square.  But  notice  that  the  covariance  vanishes, 


cov(X,X2)  =  r  xVz  dx-  (  p  xeidxj  (  f°  xV?  dx )  =0 

J — oo  \J — oo  J  \J  —  oo  / 


(9) 


Thus,  even  though  we  know  the  second  of  our  random  variables  depends  explicitly  on  the 
first,  the  standard  covariance  relationship  fails  to  help  us  understand  dependence.  While 
the  correlation  coefficient  captures  linear  dependence  between  random  variables,  we  miss 
all  the  non-linear  information. 


Copulas  address  this  problem  by  specifying  a  parametric  distribution  from  the  marginals 
of  the  random  variables.3 


To  begin  our  introduction  to  copulas,  we  require  some  definitions.  We  follow  [10]  and  [11] 
throughout. 

Definition  1  Let  S\,...,Sn  be  non-empty  subsets  in  the  extended  real  number  line,  M, 
[— oo,  oo] .  Let  H  be  a  real  function  of  n  variables  such  that  Dom  H  —  Si  x  •  •  •  x  Sn  and  for 
n-tuples  a  <  b  ( a k  <  bk  for  all  k)  let  B  —  [d.b]  =  [a\ ,b\]  x  •  •  •  x  [an,bn]  be  an  n-box  whose 
vertices  are  in  Dom  H.  Then  the  //-volume  of  B  is  given  by, 

Vtf(/?)  =  £sgn(c)//(c),  (10) 

C 

3Readers  familiar  with  Kendall’s  tau  and  Spearman’s  rho  for  testing  statistical  dependence  will  find  cou- 
plas  a  powerful  tool.  Both  Kendall’s  tau  and  Spearman’s  rho  can  be  re-expressed  as  integrals  over  copulas 
which  generate  extra  insight  into  concordance  relationships  between  random  variables.  For  further  details 
see  [10]. 
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where  the  sum  over  c  is  over  all  vertices  of  B  and  sgn(c)  denotes, 


sgn(c) 


1  if  Ck  —  cik  for  an  even  number  of  k's. 
—  1  if  ck  —  %for  an  odd  number  of  k's. 


(11) 


For  example,  if  B  is  the  3-box  [a'i  ,*2]  x  [vi  ,>’2]  x  [z\  -Zi],  then  the  //-volume  of  B  reads, 

Vh(B)  =  H(x2, yi, zi) -H(x2,y2,z\) -H{x2,y\,Z2)  -H{x\,y2,Z2) 

+  H(x2,y\,z\ )  +H(xi,y2,zi)+H(xi,yi,Z2)-H(xi,yi,zi).  (12) 

In  more  compact  notation,  we  write, 

v„(B)  =  4ff(f)  =  AS;...A5;ff(f)  (13) 

where, 

&akkH (t)  —  H(t\ ,  t2,  ■  ■  ■ ,  bk,  tk+ 1 ,  •  •  • ,  tn)  —  H(t\ ,  t2,  ■  ■  ■ ,  ak,  tk+i ,  •  •  •  ,tn).  (14) 

Definition  2  A  real  function  H  of  n  variables  is  n-increasing  if  Vh(B)  >  0  for  all  n-boxes 
whose  verticies  lie  in  Dom  H.  ■ 

We  say  that  the  function  H  is  grounded  if  H(t)  —  0  for  all  t  in  Dom  H  such  that  4.  =  ak  for 
at  least  one  k. 

Definition  3  An  n-dimensional  distribution  is  a  function  H  with  domain  in  M"  such  that  H 
is  grounded,  n-increasing  and  //(°°, . . .  ,°o)  =  1.  ■ 

If  each  Sk  has  an  greatest  element,  bk,  then  H  has  margins,  given  by 

Hk(x)  =H(bi,...,bk~i,x,bk+i,...,bn).  (15) 

We  obtain  higher  dimensional  margins  by  fixing  fewer  places. 

Definition  4  An  n-dimensional  copula  is  a  function  C  with  domain  [0, 1]"  such  that: 

•  C  is  grounded  and  n-increasing; 

•  C  has  margins  Q,  k=  1,2 which  satisfy  Ck{u)  =  u  for  all  u  in  [0,1]. 


The  above  definitions  lead  us  to  Sklar’s  Theorem: 

Theorem  1  (Sklar’s  Theorem)  Let  H  be  an  n-dimensional  distribution  function  with  mar¬ 
gins  F\ , . . . ,  Fn.  Then  there  exists  and  n-copula  C  such  that  for  all  x  in  PA , 
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H(x  i,...,xn)  =C(F\(xi),...,Fn(xn)).  (16) 

If  F\....,Fn  are  all  continuous,  then  C  is  unique,  otherwise  C  is  uniquely  determined  on 

Ran  F\  x  •  •  •  x  Ran  Fn.  Conversely,  if  C  is  an  n-copula  and  F\..  ...  Fn  are  distribution  func¬ 
tions,  then  the  function  H  defined  aboves  is  an  n-dimensional  distribution  function  with 
margins  F\ , . . . ,  Fn .  For  proof,  see  [  1 0] .  ■ 

Corollary  1  Let  H  be  an  n-dimensional  distribution  function  with  continuous  margins 
F\.....Fn  and  copula C.  Then  for  any  u  in  [0, 1]", 

C(uh  .  . . ,  un)  =  //(Ff 1  (Ml), .  . . , F~\un)).  (17) 

Again,  see  [10]  for  details.  ■ 

From  a  model  building  persepective,  copulas  allow  us  to  build  joint  distribution  functions 
from  the  marginals  of  individual  random  variables.  In  this  paper,  we  will  restrict  ourselves 
to  2-dimensional  copulas,  namely, 

H(x,y)  =  C(u,v)  =  C(F-1(x),G-1(y)),  (18) 

where  F(jt)  and  G(x)  denote  the  respective  1 -dimensional  distribution  functions  of  each 
random  variable.  We  concentrate  on  two  copula  types: 


Gaussian  Copula: 

CN(u,v;p)=<t>N{®-\u),<S>-\v);p) 


f<t>  1  (m)  /-Cl^fv) 


exp 


27t(1  —  p2)1/2  7-oo  J — oc 

where  <!>(•)  is  the  standard  normal  distribution  function. 

Gumbel  Copula  and  the  Gumbel  Survival  Copula: 


-(s2  —  2pst  +  t2 
2(1  -P2) 


dsdt, 


Cg(m,v;  a)  —  exp 


■\og(u))'l«  +  (-\og{V))'!a 


Cgs(u,  v;  a)  —u  +  v—  l-f-Qj(l  —  m,  1—  v). 


(19) 


(20) 

(21) 


We  can  easily  show  that  the  convex  sum  of  copulas  is  also  a  copula,  and  thus  we  consider 
two  model  class  for  our  foreign  exchange  analysis 

C\(u,v,0,p,a)  —  (1  —  0)Qv(«,v;p)  +  6Cg(u,v,  a) 

C2(u,v,B,p1a)  =  (1  -  e)CN(u,v;p)  +  eCGS(u,v;a)  (22) 

where  the  convexity  parameter  satisfies  0  £  [0, 1].  Using  maximum  likelihood  estimations, 
we  can  determine  which  model  best  describes  the  foreign  exchange  data. 
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Our  coupla  selection  offers  a  wide  range  of  dependency  fitting  potential  with  parsimony 
[12].  The  foreign  exchange  data  exhibits  features  contained  in  both  copula  classes.  The 
Gaussian  copula  belongs  to  a  symmetric  multivariate  elliptical  family  -  a  feature  that  for¬ 
eign  exchange  data  exhibits  at  a  gross  level  -  and  the  Gumbel  (survival)  captures  strong 
right  (left)  tail  dependence,  allowing  us  to  model  asymmetrical  fat  tails.  In  principle,  we 
could  create  a  five  dimensional  model  by  using  a  convex  sum  of  the  Gaussian,  Gumbel,  and 
Gumbel  survival  copula,  but  we  find  that  the  reduced  parsimony  does  not  give  qualitative 
better  results. 


Once  we  have  a  copula  fit  to  the  data,  we  require  a  method  for  sampling  to  generate  Monte 
Carlo  simulations.  Observe  that, 


P(V  <  v\U  =  u) 


C(u  +  Au,v)  —  C(u,v) 

A  u — >0  An 


Cu(u,v). 


(23) 


Thus  we  can  calculate  the  conditional  probability  by  taking  the  derivative  of  the  copula  with 
respect  to  one  of  its  arguments.  To  sample  from  a  copula,  we  generate  two  realizations  from 
the  uniform  distribution,  u,  and  t,  and  then  solve 


Cu(u,v)  =  t,  (24) 

for  v.  The  resulting  pair,  (w,  v),  represents  the  value  of  each  marginal  sampled  from  the  cop¬ 
ula.  Inverting  the  marginals  u  =  F(x),  and  v  =  G(y),  we  obtain  the  Monte  Carlo  realization 
of  our  random  variables. 
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4  Foreign  exchange  data  and  analysis 


We  use  ten  years  of  foreign  exchange  data  [13],  30  April  2004  to  1  May  2014,  for  USD- 
CAD,  EURCAD,  and  GBPCAD,  giving  a  total  of  2609  observations  in  each  time  series. 
Figures  1,  2,  3  show  the  historical  observations.  DND’s  largest  foreign  exchange  expo¬ 
sure  is  the  USDCAD,  followed  by  EURCAD,  and  GBPCAD  [1]  and  we  therefore  focus 
on  the  pairwise  dependency  relationship  between  the  USDCAD  and  DND’s  more  minor 
exposures.  From  a  risk  analysis  perspective,  since  the  USDCAD  dominates  DND’s  foreign 
exchange  obligations,  any  diversification  benefit  comes  primarily  from  the  USDCAD’s  re¬ 
lationship  to  the  other  two  currencies. 

In  the  QQ-plots  of  Figures  4,  5,  and  6  we  immediately  see  that  the  financial  return  data  is 
not  iidA(0, 1).  We  fit  GARCH(1,1)  models  to  our  data  using  quasi-maximum  likelihood 
methods  [8].  Using  the  Kolmogorov-Smirnov  test  on  the  filtered  returns,  we  find  that  we 
cannot  reject  the  null  hypothesis  that  the  innovations  for  each  times  belong  to  iid/V(0. 1). 
The  QQ-plots  of  Figures  7,  8,  and  9  show  the  quality  of  the  GARCH  filtered  returns  relative 
to  the  standard  normal. 

In  Figures  10,  11,  and  12  we  see  the  time  changing  stochastic  volatility  for  each  time 
series  expressed  as  an  annualized  percentage.  Notice  the  sharp  rise  in  volatility  beginning 
in  late  September  2008,  demarking  the  onset  of  the  financial  crisis.  Using  the  GARCH 
filtered  return  data,  we  compute  the  dynamics  of  the  linear  correlation  coefficient  between 
the  currency  pairs  EURCAD-USDCAD  and  GBPCAD-USDCAD;  we  display  the  results 
in  Figures  13,  and  14. 

While  the  multivariate  GARCH  model  captures  time  dependence  and  persistence,  it  does 
not  allow  asymmetrical  heavy  tailed  relationships  between  pairs  of  time  series.  To  reveal 
these  features  in  the  data,  we  use  convex  sums  of  Gaussian,  and  Gumbel  (survival)  copulas. 
We  take  each  GARCH  filtered  time  series  pair  and  convert  the  data  to  their  marginals.  We 
do  not  parametrically  model  the  marginals,  instead  relying  on  the  empirical  distribution 
functions  themselves  (recall  that  the  Kolmogorov-Smimov  test  does  not  reject  N( 0, 1)  as 
the  distribution  function  for  each  GARCH  filtered  return  time  series).  In  Figures  15,  and 
16  we  see  the  scatter  plots  in  the  marginal  space  which  immediately  show  a  Southwest  to 
Northeast  direction  -  the  time  series  pairs  exhibit  dependence,  which  we  already  discovered 
using  the  multivariate  GARCH  technique.  Fitting  the  convex  copula  sum,  we  find  that  the 
Gaussian  and  Gumbel  survival  copula  summed  together  provide  the  best  fit  to  the  data 
for  both  time  series  pairs.  We  show  in  Figures  17  and  18  the  scatter  plot  data  of  Figures 
15,  and  16  with  the  iso-contours  of  the  best  fit  copula  density  function.  The  effect  of  the 
asymmetric  tail,  captured  by  the  Gumbel  survival  copula  component,  is  readily  apparent. 
To  examine  the  robustness  of  our  result,  we  perform  a  500  sample  bootstrap  of  the  data  to 
generate  error  estimates  on  each  parameter  of  the  copula  model.  We  show  the  bootstrap 
results  in  Figure  19  and  Figure  20. 
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We  test  each  copula  model  for  time  dependence  in  each  of  its  parameters  in  turn,  holding 
the  others  fixed,  by  partitioning  the  data  according  to  a  4  x  4  grid  of  the  unit  square.  Each 
of  the  copula  parameters  0  =  (0,  p,  a),  belong  to  the  grid  through  the  sum, 

16 

©u  =  e  &tj\  (25) 

j 

where  srfj  corresponds  to  the  j-th  element  of  the  partitioned  unit  square  and  dj  is  the  16- 
tuple,  (d\,d2,  ■  ■  ■  ,6?i6),  extension  of  the  copula  parameter  considered.  For  example,  srf\  — 
[0,pi)  x  [0,gi);  s^2  —  \pi,Pi)  x  [0,<?i];  etc.  The  choice  of  the  partition  is  arbitrary  and  we 
follow  [14],  by  setting  (p,q)  G  (0.15,0.50,0.85)  to  ensure  that  we  focus  on  large  values 
within  the  partition  elements.  Thus  the  likelihood  sum  becomes, 

T 

©,-  =  argmin  J^ln(c(Mi,M2);r(zijf_1,Z2,f-i;®i)),  (26) 

t= l 

where  c(u\,  m)  is  the  copula  density  as  a  function  of  the  empirical  marginals  and, 


r(zi,f-i,z2,i-i;©0) 


(27) 


corresponds  to  the  selection  procedure  of  eq.(25).  While  this  process  can  elicit  time  depen¬ 
dence,  it  cannot  reveal  persistence  as  in  the  multivariate  GARCH  estimation  (the  partition 
method  looks  back  only  one  time  step).  The  maximum  likelihood  method  allows  us  to  esti¬ 
mate  each  dj  (16-tuple)  along  with  its  standard  error,  which  allows  us  to  test  for  statistical 
significance  in  differences  by  using  the  ^2(  1 5)  test  [14].  We  do  not  find  evidence  for  time 
dependence  in  any  of  the  copula  parameters  for  each  of  the  time  series  pairs. 

To  test  the  quality  of  the  copula  fit  to  the  data,  we  construct  a  7  x  7  contingency  table 
for  each  pairing,  comparing  the  theoretical  occupation  number  in  each  cell,  Ay,  to  the 
empirical  occupation  number  By.  The  test  statistic 


m  =  Aij)2 

i  j 


AiJ 


(28) 


is  x2  distributed  with  (k  —  l)2  degrees  of  freedom.  We  find  that  in  both  time  series  pairs, 
we  cannot  reject  the  null  hypothesis  that  the  empirical  contingency  table  is  generated  by  the 
theoretical  distribution.  We  display  the  results  in  Tables  1,  2  and  as  a  heatmap  in  Figures 
21,  22,  23,  and  24. 
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Table  1:  Seven  quantile  contingency  table  test  for  the  copula  fit  to  EURCAD  -  USDCAD. 

Quantile  counts 


Exchange  rate  pairings 

1st 

2nd 

3rd 

4th 

5th 

6th 

7th 

EURCAD-USDCAD  (empirical) 

Quantile  counts 

7th 

16 

21 

26 

38 

46 

74 

151 

6th 

18 

29 

42 

48 

74 

80 

82 

5th 

30 

40 

55 

63 

81 

60 

44 

4th 

38 

51 

68 

82 

57 

51 

25 

3rd 

48 

77 

57 

52 

54 

52 

33 

2nd 

83 

77 

72 

54 

35 

32 

20 

1st  139  78 

EURCAD-USDCAD  (predicted) 

53 

35 

26 

24 

17 

Quantile  counts 

7th 

10 

20 

29 

40 

54 

79 

139 

6th 

20 

33 

41 

52 

66 

81 

79 

5th 

29 

41 

50 

62 

70 

66 

54 

4th 

39 

50 

62 

69 

62 

52 

40 

3rd 

52 

65 

73 

62 

50 

41 

29 

2nd 

76 

88 

65 

50 

41 

33 

20 

1st 

146 

76 

52 

39 

29 

20 

10 

Z2(36) 


43.4  (c.v.  51.0) 
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Table  2:  Seven  quantile  contingency  table  test  for  the  copula  fit  to  GBPCAD  -  USDCAD. 

Quantile  counts 


Exchange  rate  pairings 

1st 

2nd 

3rd 

4th 

5th 

6th 

7th 

GBPCAD-USDCAD  (empirical) 

Quantile  counts 

7th 

9 

16 

20 

33 

50 

89 

155 

6th 

13 

28 

44 

56 

63 

80 

89 

5th 

19 

35 

49 

67 

80 

59 

63 

4th 

39 

52 

63 

62 

82 

51 

24 

3rd 

48 

71 

65 

69 

54 

40 

26 

2nd 

84 

88 

75 

51 

28 

36 

11 

1st  160  83 

GBPCAD-USDCAD  (predicted) 

57 

34 

16 

18 

4 

Quantile  counts 

7th 

5 

14 

24 

37 

55 

85 

154 

6th 

14 

28 

41 

54 

69 

83 

85 

5th 

23 

40 

53 

63 

70 

69 

55 

4th 

35 

53 

64 

67 

63 

54 

37 

3rd 

52 

69 

72 

64 

53 

41 

24 

2nd 

82 

87 

69 

53 

40 

28 

14 

1st 

163 

82 

52 

35 

23 

14 

5 

Z2(36) 


36.1  (c.v.  51.0) 
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Figure  1:  USDCAD  performance  April  30,  2004  -  May  1,  2014. 
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Figure  2:  EURCAD  performance  April  30,  2004  -  May  1,  2014. 


14 


DRDC-RDDC-2014-R62 


Figure  3:  GBPCAD  performance  April  30,  2004  -  May  1,  2014. 
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Figure  4:  USDCAD  QQ-plot  returns  vs  standard  normal  April  30,  2004  -  May  1,  2014. 
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Figure  5:  EURCAD  QQ-plot  returns  vs  standard  normal  April  30,  2004  -  May  1,  2014. 
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Figure  6:  GBPCAD  QQ-plot  returns  vs  standard  normal,  April  30,  2004  -  May  1,  2014. 


18 


DRDC-RDDC-2014-R62 


Figure  7:  USDCAD  QQ-plot  GARCH  filtered  returns  vs  standard  normal  April  30,  2004  - 
May  1,  2014. 
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Figure  8:  EURCAD  QQ-plot  GARCH  filtered  returns  vs  standard  normal  April  30,  2004  - 
May  1,  2014. 
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Figure  9:  GBPCAD  QQ-plot  GARCH  filtered  returns  vs  standard  normal,  April  30,  2004 
-May  1,  2014. 
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Figure  10:  USDCAD  GARCH  modelled  annualized  volatility  April  30,  2004  -  May  1, 
2014. 
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Figure  11:  EURCAD  GARCH  modelled  annualized  volatility  April  30,  2004  -  May  1, 
2014. 
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Figure  12:  GBPCAD  GARCH  modelled  annualized  volatility,  April  30,  2004  -  May  1, 
2014. 
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Figure  13:  The  correlation  coefficient  for  the  EURCAD-USDCAD  pair  April  30,  2004  - 
May  1,  2014. 
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Figure  14:  The  correlation  coefficient  for  the  GBPCAD-USDCAD  pair  April  30,  2004  - 
May  1,  2014. 
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Figure  15:  EURCAD-USDCAD  marginals,  April  30,  2004  -  May  1,  2014. 
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Figure  16:  GBPCAD-USDCAD  marginals,  April  30,  2004  -  May  1,  2014. 
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Figure  17:  EURCAD-USDCAD  Gaussian  Gumbel  survival  copula  density  fit,  April  30, 
2004 -May  1,  2014. 
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Figure  18:  GBPCAD-USDCAD  Gaussian  Gumbel  survival  copula  density  fit,  April  30, 
2004 -May  1,  2014. 
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Figure  19:  EURCAD-USDCAD  bootstrapped  (500  samples)  copula  parameter  error  distri¬ 
butions. 


DRDC-RDDC-2014-R62 


31 


Figure  20:  GBPCAD-USDCAD  bootstrapped  (500  samples)  copula  parameter  error  distri¬ 
butions. 
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Figure  21:  Heatmap  depiction  of  the  EURCAD-USDCAD  empirical  contingency  table. 
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Figure  22:  Heatmap  depiction  of  the  EURCAD-USDCAD  copula  predicted  contingency 
table. 
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Figure  23:  Heatmap  depiction  of  the  GBPCAD-USDCAD  empirical  contingency  table. 
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Figure  24:  Heatmap  depiction  of  the  GBPCAD-USDCAD  copula  predicted  contingency 
table. 
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Table  3:  100  trading  day  VaR  levels  from  May  1,  2014  for  a  portfolio  short  3  million  USD 
and  1  million  GBR  Results  based  on  Monte  Carlo  simulation  (sample  size:  3  x  104). 


Model 

VaR  io%  ($M  CAD) 

VaR5%  ($M  CAD) 

VaRi%  ($M  CAD) 

Copula-GARCH 

0.24 

0.32 

0.48 

Multivariate-GARCH 

0.27 

0.34 

0.52 

Independence 

0.21 

0.27 

0.41 

5  VaR  analysis 


We  can  use  the  copula  and  multivariate  GARCH  models  to  build  VaR  estimates  through 
Monte  Carlo  techniques.  The  multivariate  GARCH  approach  ensures  that  the  resulting 
simulations  contain  the  modelled  correlation’s  time  dependence  and  persistence  between 
exchange  rate  pairs.  The  copula  GARCH  method  models  the  joint  distribution  function  as 
a  fixed  copula,  ensuring  VaR  levels  consistent  with  dependency  relationships. 

To  demonstrate  the  applicability  of  the  VaR  estimate,  consider  a  hypothetical  project  that  is 
short  3  million  USD  and  short  1  million  GBP,  both  due  100  trading  days  after  May  1,  2014. 
In  table  3  we  see  the  foreign  exchange  VaR  levels  relative  to  the  May  1,  2014  spot  purchase 
of  5.14  million  CAD.  We  see  that  both  dependency  models  give  similar  results  while  the 
independence  assumption  underestimates  the  risk  in  the  portfolio:  at  the  1%  VaR  level,  the 
independence  assumption  yields  a  result  25%  lower  than  the  dependency  models.  In  figure 
25  and  26  we  show  1000  simulated  joint  exchange  rate  realizations  for  both  dependency 
models. 
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Figure  25:  GBPCAD-USDCAD  simulation:  copula  scatter  plot  100  trading  day  horizon 
from  May  1,  2014  (1000  points). 
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Figure  26:  GBPCAD-USDCAD  simulation:  multivariate  GARCH  scatter  plot  100  trading 
day  horizon  from  May  1,  2014  (1000  points). 
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6  Discussion 


We  present  two  methods  -  the  multivariate  GARCH,  and  the  copula-GARCH  -  for  esti¬ 
mating  the  foreign  exchange  return  process.  Using  both  methods,  we  can  compute  the  VaR 
for  a  portfolio  of  foreign  exchange  obligations,  properly  including  the  dependence  between 
currency  pairs.  We  can  construct  Monte  Carlo  samples  from  each  model  and  reverse  the 
GARCH  process  to  give  simulated  returns.  Both  models  offer  advantages  -  the  multivari¬ 
ate  GARCH  model  captures  the  persistence  of  shocks  to  the  correlation  coefficient,  while 
the  copula-GARCH  model  captures  asymmetric  tail  dependence.  Using  both  models  will 
help  DND  managers  understand  multiple  foreign  currency  exposures  risk  within  project 
budgets. 

We  find  that  the  pairings  EURCAD-USDCAD  and  GBPCAD-USDCAD  exhibit  strong 
dependence  in  both  modelling  approaches.  When  the  Canadian  dollar  does  well  or  poor 
against  one  of  the  currencies,  the  Canadian  dollar  tends  to  follow  the  same  performance 
against  the  other  currency.  Thus,  DND’s  multiple  currency  exposures  provides  some,  but 
limited,  diversification  benefit  -  especially  during  tail  events. 

Senior  decision  makers  can  use  these  models  to  gauge  the  foreign  exchange  risk  within 
project  budgets.  We  see  that  the  exposure  to  multiple  currencies  does  not  eliminate  foreign 
exchange  risk  and  we  cannot  simply  convolve  distributions  from  univariate  time  series  for 
each  currency  to  estimate  portfolio  VaR.  The  dependency  between  the  exchange  rates  plays 
an  important  role  in  determining  the  overall  VaR  level  of  the  portfolio. 

While  VaR  offers  a  risk  monitoring  mechanism  for  DND,  it  does  not  in  and  of  itself  offer 
risk  mitigation.  DND  must  decide  on  the  risks  it  is  comfortable  holding.  As  a  cautionary 
note,  before  embarking  on  any  risk  mitigation  strategy  in  foreign  exchange  -  contingency 
reserve  or  derivative  based  -  DND  must  understand  the  source  of  possible  gains.  In  fi¬ 
nancial  firm  evaluation,  creating  or  buying  insurance  against  unfavourable  outcomes  does 
not  affect  the  value  of  a  firm.  Insurance  repackages  the  firm’s  cash  flows,  and  cash  flow 
repackaging  cannot  increase  firm  value.4  If  insurance  does  provide  value,  it  must  come 
from  advantages  that  the  insurance  separately  creates.  For  example,  if  management  can 
re-focus  its  efforts  in  areas  where  it  has  a  comparative  advantage  by  entering  an  insurance 
contract,  then  the  repackaging  can  have  a  positive  effect.  In  DND’s  case,  financial  risk 
mitigation  of  any  kind  must  be  tied  to  the  real  source  of  potential  gains  -  superior  project 
management. 

CORA  has  afforded  DND  three  powerful  methods  for  modelling  and  understanding  foreign 
exchange  risk  across  department  budgets:  mulitvariate  GARCH  models,  copula-GARCH 
models,  and  models  based  on  the  prices  of  vanilla  foreign  exchange  options  (in  particular, 
the  risk  neutral  probability  distribution  from  the  Heston  model).  Each  modelling  method 

4This  is  a  statement  of  the  famous  Modigliani-Miller  Theorem  -  a  no  arbitrage  argument  for  firm  pricing. 
See  [15]  for  details. 
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offers  insight  into  the  exchange  rate  market  risk  that  DND  assumes  in  foreign  transactions. 
They  capture  realistic  features  of  market  data  including  heavy  tails,  asymmetric  density 
functions,  and  leptokurtosis.  The  inclusion  of  regular  VaR  updates  in  risk  reporting  over 
time  horizons  sensitive  to  major  project  deliveries  will  help  senior  DND  decision  makers 
better  understand  their  foreign  currency  liabilities  in  military  operations. 
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